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ABSTRACT 

We present multi-wavelengths observations and a radiative transfer model of a newly discov- 
ered massive circumstellar disk of gas and dust which is one of the largest disks known today. 
Seen almost edge-on, the disk is resolved in high-resolution near-infrared (NIR) images and ap- 
pears as a dark lane of high opacity intersecting a bipolar reflection nebula. Based on molecular 



^ Based on observations made at the Calar Alto Observatory. This worlc is based in part on data collected at Subaru 
Telescope, which is operated by the National Astronomical Observatory of Japan, and on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with 

NASA. 



line observations we estimate the distance to the object to be 3.5 kpc. This leads to a size for 
the dark lane of ~10500 AU but due to shadowing effects the true disk size could be smaller. 
In Spitzer/IRAC 3.6 fim images the elongated shape of the bipolar reflection nebula is still pre- 
served and the bulk of the flux seems to come from disk regions that can be detected due to the 
slight inclination of the disk. At longer IRAC wavelengths, the flux is mainly coming from the 
central regions penetrating directly through the dust lane. Interfcrometric observations of the 
dust continuum emission at millimeter wavelengths with the SMA confirm this finding as the 
peak of the unresolved mm-emission coincides perfectly with the peak of the Spitzer/IRAC 5.8 
fm\ flux and the center of the dark lane seen in the NIR images. Simultaneously acquired CO 
data reveal a molecular outflow along the northern part of the reflection nebula which seems to 
be the outflow cavity. An elongated gaseous disk component is also detected and shows signs of 
rotation. The emission is perpendicular to the molecular outflow and thus parallel to but even 
more extended than the dark lane in the NIR images. Based on the dust continuum and the CO 
observations we estimate a disk mass of up to a few solar masses depending on the underlying 
assumptions. Whether the disk-like structure is an actual accretion disk or rather a larger-scale 
flattened envelope or pseudodisk is difficult to discriminate with the current dataset. The exis- 
tence of HCO+/H^'^CO"'" emission proofs the presence of dense gas in the disk and the molecules' 
abundances are similar to those found in other circumstellar disks. We furthermore detected C2H 
towards the objects and discuss this finding in the context of star formation. Finally, we have 
performed radiative transfer modeling of the K band scattered light image varying a disk plus 
outfiow 2D density profile and the stellar properties. The model approximately reproduces extent 
and location of the dark lane, and the basic appearance of the outflow. We discuss our findings 
in the context of circumstellar disks across all mass regimes and conclude that our discovery 
is an ideal laboratory to study the early phases in the evolution of massive circumstellar disks 
surrounding young stellar objects. 

Subject headings: stars: pre-main sequence, stars: formation, planetary systems: protoplanetary disks, 
infrared: stars 



1. Introduction 

In recent years, significant progress has been 
made in our understanding of Young Stellar Ob- 
jects (YSOs) and the formation of stars. Still, 
there remain unanswered questions and problems 
that have to be solved observationally and the- 
oretically. In particular, while the formation of 
low-mass and intermediate mass stars by accre- 
tion of material via a circumstellar disk is well 
established (e.g., Watson et al.||2007 ), for massive 
stars (M > 8 Mq) not only the theory of the star 
formation process (e.g., iStahler et al.||2000|) but 



also the observations are challenging (e.g., Beuther 
et al.||2007| [Zinnecker fc Yorkel|2007| a nd yield 



partly contradicting results (e.g., Chini et al._2004l 



Sako et al.|2005l|Steinacker et al.|2006[|Niirnberg'CT| 



et al.|2007 ) . Conceptually the problem arises from 
the fact that at least in the spherically symmetric 
case the strong radiation pressure of a hydrogen 
burning young massive star with a mass > 8 M© 
would prevent any further gas infall and would 
thus halt the formation of a more massive object. 



Two competing theories overcoming this problem 
were intensively discussed in the recent past. They 
explain the formation of massive stars either by 
the so-called " Competitive Accretion and Merging 
Model" (iBonnell & Bate|'2006) or by fragmenta- 



tion during a large scale collapse in turbulent cores 
(iMcKee & Tan||2003|). The latter model resem- 



bles the formation process of low-mass objects and 
consequently relies on circumstellar disks around 
young massive stars to build up the final objects 
(see, also, Krumholz et al.pOOQ ). 

Today, there is ample evidence for the ex- 
istence of "disk- like" structures around mas- 
sive young stellar objects with masses/sizes be- 
ing higher/bigger compared to typical low-mass 
circumstellar disks (e.g., iCesaroni et al.| |2007 



et al 



2004| Patel et al.| 2005 Jiang et al. 2005 



Schreyer et al.| [20061 |Chini et al.||2004[ [Beltran 



2008]). While for early B-stars these structures 
appear to be indeed rather disk-like, for 0-stars 
they appear to have more of a toroidal shape. 



In most cases large, flattened structures were de- 
tected using high spatial resolution molecular line 
data in the (sub-)nim regime. Typically, at shorter 
wavelengths, the dense envelope surrounding the 
massive protostar causes too much extinction and 
prevents a direct detection of the source in the 
mid- or near-infrared. Indications for rotation 
in the flattened gaseous structures supported the 
hypothesis that these objects are massive circum- 
stellar disks. However, the question whether these 
objects are indeed the high-mass counterpart to 
the well-established low-mass protoplanetary ac- 
cretion disks directly feeding the central source or 
whether they are (remnants of) the rotating and 
inf ail ing molecular envel opes (e.g., Beuther et al. 

( |GalhfcShu||1993| [Allen 



MIR images, Chini et al. (20061 reported the de- 



2009) or "pseudodisksl 



et al. 2003) is unclear. In particular, it appears 



likely that at least some of the observed massive 
gaseous "toroids" around stars of 10-20 Mq are 
unstable and have comparatively short lifetimes 



(Cesaroni et al. 2007) while low-mass disks typi- 



cally persist for several Myr. 

In addition to (sub-)nim data, only for very few 
massive objects mid- or near-infrared (MIR / NIR) 
data exist that support the existence of an ex- 



tended dusty disk component. Chini et al. (2004) 



presented an NIR image of an assumed accretion 
disk surrounding a presumably massive protostar 
in M17, but their interpretation was challenged 



by Sako et al. ( 2005 1 . A detailed radiative trans- 



fer modeling of the iiT-band image of this source 
was performed in Steinacker et al. (2006) show- 
ing that the underlying structure has a disk shape 
with a mass around a few M© . Using NIR polari- 



metric images Jiang et al. ( 2005 2008 ) found elon 



gated polarized signatures around massive young 
stellar objects suggesting the presence of polar- 
ization disks. Sridharan et al. (2005) claimed 



the detection of a silhouette edge-on disk around 
IRAS 20126-F4104 based on NIR images in the K, 
L', and M' filters. More recently, however, De 



Buizer (2007) derived from MIR images at 12.5 



and 18.3 /zm a much more complex structure of 
the inner regions of the same source and suggested 
that the data could also be interpreted having a 
tight young stellar cluster in the center. Based on 



^A " pseudodisk" is a flattened disequilibrium structure cre- 
ated during the collapse of a magnetized cloud core by 
strong magnetic pinching for ces deflecting the i nfalling gas 
toward the equatorial plane (Galli & Shu|1993L 



tection of a remnant disk around a young massive 



star, and Nielbock et al. (2007) reported on evi 



dence for a disk candidate surrounding a hyper- 
compact HII region in M17 based on NIR data. 
These authors also provided a radiative transfer 
model for the disk seen in the K^-band image. Fi- 



nally, very recently Okamoto et al. (2009) pub- 



lished the detection of a flared circumstellar disk 
around the 10 Mq Herbig Be star HD200776 also 
using MIR imaging data. 

Here, we report the direct detection of a large 
and massive circumstellar disk of dust and gas sur- 
rounding a young stellar source. The object is 
named CAHA J23056-I-6016 as it was first iden- 
tified in NIR images taken at the Calar Alto 
observatory (CAHA). Currently, we can not say 
whether the disk is an actual accretion disk or 
rather a "pseudodisk" or a flattened envelope and 
we use the more general term " circumstellar disk" 
or simply "disk" throughout the paper to avoid 
any possible misinterpretation. The morphology 
shows a bipolar nebula intersected by a dark lane. 
Multi-wavelengths data including NIR/MIR im- 
ages as well as dust continuum and molecular line 
data are used to constrain the physical properties 
of this object. We derive a distance of 3.5 kpc 
to the source, find a radius for the disk of several 
thousand AU and a mass, depending on the under- 
lying assumptions, of up to several solar masses. 
The K—hand and the mm image are analyzed with 
a disk/outfiow model using radiative transfer cal- 
culations to derive scattered light images. 

In section 2 we describe the observations and 
data reduction and present the results of the data 
in section 3. Results from the radiative trans- 
fer model are presented in section 4. A discus- 
sion about the nature and properties of CAHA 
J23056+6016 in the context of star formation is 
given in section 5. Finally, we summarize our find- 
ings and conclude in section 6. 

2. Observations and data reduction 
2.1. Omega2000 NIR wide-field images 

CAHA J23056-I-6016 was detected in a deep 
NIR survey in October 2003 with the wide-field 
NIR camera Omega2000 on the 3.5 meter tele- 



scope in Calar Alto (SpainFJ The camera is 
equipped with a 2048x2048 pixel HAWAII-2 de- 
tector, has a pixel scale of 0.45"/pixel and an ef- 
fective field-of-view of ^15'. Images were taken in 
the three NIR filters J , H , and Kg with central 
wavelengths at 1.21, 1.65, and 2.15 /xm, respec- 
tively. The observations followed a pre-defined 
dither pattern of 20 positions around a defined 
image center. On each position 20 frames with 
an integration time of 3 seconds were taken, sum- 
ming up to 1 minute integration time per dither 
position. After co-addition of all images the to- 
tal integration time amounted to 20 minutes per 
filter. 

Flat fielding, bad pixel correction and co- 
addition of the images were done with standard 
IRAFg routines. Furthermore we used the driz- 
zle routine to reduce the pixel scale of the fi- 
nal images to half of its original value, i.e., to 
0.225"/pixel. Astrometry was applied by compar- 
ing the positions of detected 2MASS point sources 
to those listed in the 2MASS Point Source Cat- 
alog (PSC) ( |Cutri et al!]|2003[ ). On average the 
precision of the astrometry is better than 0.25" 
with respect to the known 2MASS positions. 

The photometric calibration of CAHA J23056+6016 
and its nebula (see section 3.1) was done using 
six 2MASS sources in the immediate vicinity of 
CAHA J23056-I-6016 as reference stars (for the 
J filter image seven reference sources were avail- 
able). The IDL routine ATV.pro was used to 
carry out aperture photometry and derive count 
rates of the sources in the Omega2000 images. 
The 2MASS magnitudes of the objects were then 
converted to mJy using the Magnitude-to-Flux 
Density converter provided by the Spitzer Science 
CenteiQ Plotting flux density vs. count rate for 
all reference object and fitting an error-weighted 
line though the data allowed us to derive the re- 
lation between counts/pixel and mJy/per pixel. 
For each filter this relation was used to calibrate 
our Omega2000 images. As a cross-check we per- 
formed aperture photometry on the calibrated im- 
ages to re-derive the fiux density for the reference 
sources. For each filter the results were within the 
error bars from the 2MASS photometry giving us 



confidence in our photometric calibration. 

2.2. Subaru IRCS high-resolution NIR im- 
ages 

In June 2008 high-resolution H- and iiT'-band 
images were obtained with the Infrared Camera 
and Spectrograph (IRCS) at the Subaru telescope. 
IRCS is equipped with a 1024x1024 pixel Al- 
ladin3 detector providing a pixel scale of ~0.027 
mas/pixel. In each filter 7 sets of 9 dither posi- 
tions were observed. On each dither position 5 
exposures of 5 seconds were taken resulting in a 
total integration time of 1575 sec (=26.25 min) 
per filter. The data processing and photometric 
calibration was done in the same manner as we 
did for the wide-field images. We reached a 5-ct 
limiting magnitude of ~22.7 mag in H and ^-^23.1 
mag in K in the 2MASS photometric system. 

2.3. Spitzer IRAC images 

Spitzer IRAC data were downloaded from the 
Spitzer data archive (AOR-Key: 18917120). The 
data were taken in February 2007 and serendipi- 
tously covered CAHA J23056-I-6016. As the ob- 
serving parameters were not optimized for CAHA 
J23056-I-6016 but rather for bright, nearby HII 
regions, the integration time per pixel was only 
1.2 seconds. Hence, for CAHA J23056-I-6016 the 
data are of modest signal-to-noise and only data 
from channel 1, 2 and 3 (at 3.6, 4.5 and 5.8 /im, 
respectively) are used in our analyses. We used 
M0PEX[J to carry out the basic data reduction 
steps (outlier rejection, mosaicking) and create the 
final images. In this process, the pixel scale was 
changed to 0.6"/pixel in all channels instead of the 
default value of 1.2"/pixel. To convert the flux 
densities from the initially given unit of MJy/sr 
to /xJy/pixel (similar to the NIR data) we adopted 
the value provided by the Spitzer Science Centeir] 
to the new pixel scale: 



1 MJy/sr « 8.46161 ^Jy/pixel 



(1) 



^http://ssc. spitzer. caltoch. edu/postbcd/download- 

mopex.html 
^http://ssc. spitzer. caltech.edu/archanaly/quickphot. html 



2http://w3. caha.es/CAHA/Instruments/O2000/index2. html 

•^ http://iraf.noao.edu/ 

*http://ssc. spitzer. caltech.edu/tools/magtojy/ 



2.4. Millimeter wavelength observations 

2.4.1. SMT 

Following the initial discovery of CAHA J23056+6016 
in 2003 we obtained first radio data at the 10-m 
Sub-Millimeter Telescope (SMT) on Mount Gra- 
ham in December 2004. Molecular line measure- 
ments of CO(3-2) (beam size ~22") were obtained 
at 345 GHz. We used the facility receivers and re- 
duced all spectra with the CLASS software^ The 
atmospheric transmission was T225 ghz ~ 0.2 and 
W75N was used as calibration source. The esti- 
mated pointing error is 3" and the radiometric 
accuracy is 15%. 

2.4.2. IRAM 30m telescope 

In June 2006 we were granted DDT time at 
the IRAM 30-m telescope on Pico Veleta (Spain) 
to observe CAHA J23056-I-6016 in mm dust con- 
tinuum with MAMBOl in on-off mode. The 
MAMBOl bolometer works at 1.2 mm, has 37 
pixels and a pixel spacing of 20" . The half-power- 
beam-width (HPBW) is -10.5". The total on- 
source integration time amounted to 40 minutes. 
Cep A and NGC7538 served as calibration sources. 
The data reduction was done using standard rou- 
tines from the MOPSIC software package. 

Furthermore, the IRAM 30 m telescope was 
used to obtain spectra of CAHA J23056+6016 
in the molecular lines HCO+(1-0) at 89.189 GHz 
and its isotopologue }1^^CO+(1-0) at 86.754 GHz. 
These observations took place between August 20 
- 24, 2009. The EMIR receiver E90 was utihzed in 
combination with the VESPA backend. Position 
switching was employed with on-source times of 
120 s and 840 s for HCO+(1-0) and Hi3C0+(l- 
0), respectively. The spectral resolution was 40 
kHz per channel, which corresponds to roughly 
0.135 km/s. The data were first calibrated to 
the antenna temperature scale using the common 
chopper-wheel method and then converted into 
main beam brightness temperatures using the for- 
ward and main beam efficiencies provided by the 
observatorjFJ The achieved rms noise levels (in 
Tmb) are 0.043 K for HCO+(1~0) and 0.018 K for 
H"CO+(1-0), respectively. 



The C2H(l-0) transitions were observed in July 
2009 at the IRAM 30 m telescope with the new re- 
ceiver system EMIR. All hyperfine structure lines 



near 87.4 GHz were observed (Bel & Leroy 1998) 
with a spectral resolution of ~0.13kms~^. Sys- 
tem temperatures at 3 mm wavelengths were 88 K. 
The la rms of the final smoothed spectrum with 
0.27kms~l spectral resolution is 25 mK. 

2.4.3. SMA 



CAHA J23056-f6016 was observed with the 
SubmiUimeter Array (SMAJj^on October 7th 2007 
with 8 antennas at 1.3 mm in the compact con- 
figuration with projected baselines up to 53 kA. 
The phase center was 23''05™37.4" in RA, and 
60°15'45".6 in DEC (J2000.0) with a tuning fre- 
quency of 230.537970 GHz in the upper sideband 
and chunk 14 (uigi. = — 49.6kms~"'^). The spectral 
resolution was 0.406 kHz/channel corresponding 
to ~ 0.55 km s~^ velocity resolution. The weather 
was intermediate with zenith opacities r(225GHz) 
around 0.15 measured by the National Radio As- 
tronomy Observatory (NRAO) tipping radiometer 
operated by the Caltech SubmiUimeter Observa- 
tory (CSO). Passband and flux calibration were 
derived from 3C111 and Ganymede observations, 
and the flux density scale is estimated to be accu- 
rate within 20%. Phase and amplitude were cal- 
ibrated with regularly interleaved observations of 
the quasar 0102-1-584. We applied different weight- 
ings for the continuum and line data resulting in 
synthesized beams of 4.3" x 3.7" for the mm con- 
tinuum emission and 3.8" x 3.3" and 4.0" x 3.5" for 
the Ci^0(2-1) and i2C0(2-l) emission, respec- 
tively. The rms noise of the line and continuum 
images are 1.4mJybeam~^ and 150mJybeam~^ 
per 0.6kms~^ wide channel. The initial flagging 
and calibration was done with the IDL superset 
MIR originally developed for the Owens Valley Ra- 
dio Observatory ( Scoville et al.|1993 ) and adapted 
for the SMAp°| The imaging and data analysis 
were conducted in MIRIAD ( |Sault et alT]|1995 |. 



^http://www.iram.fr/IRAMFR/GILDAS 

*http://www. irain.es/IRAMES/mainWiki/Iram30mEfHciencies 



^The SubmiUimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics, and is 
funded by the Smithsonian Institution and the Academia 
Sinica. 
^"The MIR cookbook by Charlie Qi can be found at 
|http://cfa- www. harvard. edu/~cqi/inircook.html| 



3. Results 

3.1. NIR images 

3.1.1. The large scale environment 

Figure [l] shows a color composite of the 
OMEGA2000 NIR images with CAHA J23056+6016 
lying at the image center (RA (J2000): 23h 05m 
37.52s, DEC (J2000): +60° 15' 45.8"). A bipolar, 
cone-shaped nebula and an intersecting dark lane 
are clearly visible. In the extended surroundings 
of CAHA J23056+6016 there is evidence for ongo- 
ing star formation: The large HII complex in the 
lower right corner of the image is IC 1470 (alias 
IRAS 23030-^5958 or S156) with a distance of 
~ 3.6' from CAHA J23056+6016, and the bright 
star in the lower left corner is WRAM17 (distance 
~ 2.7') which is assumed to be the exciting source 
of BFS17, a compact HII region in its immedi- 
ate vicinity (Russeil et al. 20071. Furthermore, 



IRAS 23033-1-5951, another compact high-mass 
star-forming region, is lying ~ 7.7' south of CAHA 
J23056-H6016. 



3.1.2. Source morphology 
disk 



an almost edge-on 



Figure [2] focuses directly on the detailed struc- 
ture of CAHA J23056-I-6016 in the NIR. While 
the J filter image was taken at Calar Alto, the H 
and K' filter images were obtained with Subaru 
and superior spatial resolution. The contours give 
an idea of the different observed fiux levels and 
the achieved S/N in each filter. Assuming that 
the flux is mostly scattered light from a dusty sur- 
face (see below) the fact that the nebula appears 
brighter at longer NIR wavelengths is surprising 
as the scattering efficiency increases for shorter 
wavelengths for typical ISM dust grains. However, 
line-of-sight extinction towards the source causes 
higher absorption at short wavelengths and could 
explain the difference in the apparent brightness 
of the nebula. 

The high-resolution images in the middle and 
right-hand panel in Figure [2] clearly resolve the 
bipolar nebula. The nebula is rotated in the plane 
of the sky by roughly ~10° from north to the west 
and stretches over ~12" in the north-south and 
>6" in the east-west direction. The northern and 
southern cones of the nebula are asymmetric in 
size and brightness, with the northern part being 



more prominent in size and luminosity. Assuming 
an intrinsic symmetry of the structure this can 
be explained if the northern part is tilted towards 
the observer compared to the southern part (see 
also sections 3.2 and 3.3.3) and if forward scatter- 



ing from dust grains is the main source of bright- 
ness of the nebula in the NIR. Forward scattering 
and the suggested orientation would also explain 
the rather smooth brightness distribution and the 
position of the brightest regions in the nebula. 
For the scattering to work, an internal luminos- 
ity source must be present illuminating its dusty, 
cone shaped surroundings. The central source, 
however, is not detected in the NIR images di- 
rectly but is presumably hidden behind the flat- 
tened structure of high extinction intersecting the 
northern and southern part of the nebula. As can 
be seen in Figure [2] the degree of extinction of this 
intersecting lane depends on the wavelengths as 
would be expected if dust was the main source of 
extinction: Even if the J-band image has a lower 
S/N than the other images the gap between the 
northern part of the nebula and the southern part 
appears much wider. Going to longer wavelengths, 
the apparent size of the gap decreases in H, and 
with the chosen contour levels the gap is no longer 
apparent in the K' image. However, as can be 
nicely seen in the K' band, the contours suggest 
a layered structure of extinction having its maxi- 
mum in the middle and then slowly decreasing to- 
wards the north and the south. Putting all this to- 
gether, CAHA J23056-I-6016 shows a striking sim- 
ilarity to edge-on circumstellar disk systems found 
around young low-mass stars with a characteristic 
bipolar nebula and an intersecting dark dust lane 



(e.g., Padgett et al. 1999 Brandner et al. 20001 



However, as mentioned above, due to the asymme- 
try and rotation of the nebula the disk of CAHA 
J23056-f6016 appears to be slightly inchned and 
tihed. 

The apparent size of the dark lane at the base 
of the northern part of the nebula is ~3". This 
translates into ^10500 AU for an assumed dis- 
tance of 3.5 kpc (see section [3.3.1 for details on 
the derivation of the distance). However, the dark 
lane is caused by the extinction of scattered light 
in the disk and what we observe could actually be 
not the disk itself but rather its shadow onto the 



surrounding nebula (Pontoppidan & Dullemond 



2005). Hence, in this case, we can not derive the 



disk size directly from the images and as described 
in section 4.1 the actual disk size may indeed be a 
bit smaller. 

The nebula in Figure [2] has a physical size of 
~ 42000 X 21000 AU or ~ 0.2 x 0.1 pc and ap- 
pears to be more than one order of magnitude 
larger than those detected for edge-on systems in 



nearby star-forming regions (e.g., Padgett et al. 
T9991 [Brandner et al.||2000[ ). 



3.2. IRAC images - peaking through the 
dust lane 

Going to slightly longer wavelengths, Fig- 
ure Is] shows Spitzer/IRAC images of CAHA 
J23056-I-6016. As mentioned in section 2.3. no 
optimized pointed observations were carried out 
but CAHA J23056-F6016 was rather serendipi- 
tously imaged during a larger mapping program. 
Hence, the S/N is limited, especially in the images 
at 3.6 and 4.5 fim (see figure caption for details). 

At 3.6 /im the morphology of the emission still 
slightly resembles the overall morphology in the 
NIR. The emission is elongated in north-south di- 
rection and appears to be rotated slightly west- 
wards. At this moment, we do not know whether 
the observed flux at this wavelength is mainly 
coining from PAH emission or whether it is caused 
by scattering. At 4.5 /im the emission becomes al- 
ready more compact and the rotation of the elon- 
gated structure is gone. As mentioned above the 
emission in the third channel at 5.8 /im is rather 
intense and, similar to channel 2, the morphology 
does no longer resemble the NIR images. Two 
things are noteworthy regarding the IRAC images: 
(1) Comparing the peak of the emission in each 
channel one finds that it apparently changes its 
position. Coming from 3.6 /im the emission peak 
moves slightly southwards at 4.5 fim and then fur- 
thermore to the west at 5.8 /im. The overall shift 
between IRAC 1 and 3 is shown in Figure[2]and ap- 
pears to be > 0.8". This is significantly more than 
the expected astrometric accuracy of post-BCEJ^ 
IRAC data which is better than 0.3" (see, Spitzer 
webpaga^. Without emphasizing to much the 
derived numbers we believe that at least the trend 
is real. (2) Compared to the maximum emission 
in the NIR which was confined to the nebula, the 



^^post Basic Calibrated Data 

-"^^httpiZ/ssc. spitzer. caltech.edu/irac/products/overview. html 



maximum emission in the IRAC channel 2 and 3 
is centered directly on the dark lane. 

Building upon the dust disk hypothesis from 
the NIR images these findings can be interpreted 
as follows. At longer IRAC wavelengths the disk 
is no longer opaque and one starts to detect the 
inner regions of the object potentially combining 
emission from the inner disk regions and the cen- 
tral source. 

The offset in the emission peak at 3.6 /im com- 
pared to 5.8 /xm can be explained if one assumes, 
as laid out in section 3.1.2., that the disk is not 
seen edge-on but slightly inclined and rotated in 
the plane of the sky. The maximum flux at 3.6 
/xm would then not come from the central source 
because the opacity of the disk is still too high to 
directly " see" through the disk at this wavelength. 
Rather, the emission presumably traces the heated 
inner disk rim and those parts of the disk surface 
that are tilted towards the observer and are sub- 
ject to less extinction. In Figure [2] it can be seen 
that the emission peak at 3.6 /im is located at a 
position of higher luminosity and thus less extinc- 
tion than the emission peak at 5.8 /im. Figure [4] 
shows a color-composite of H, K' and IRAC 4.5 
/im where the different regions of maximum emis- 
sion become nicely evident. 

3.3. SMT, IRAM and SMA (sub-)mm 
data 

3.3.1. Distance estimate 

To derive further information about CAHA 
J23056-I-6016 and to probe the properties of dust 
and gas in its vicinity, additional observations at 
even longer wavelengths were carried. In Figure [5] 
we show the result for the CO (3-2) molecular 
line observations at the SMT with the 345 GHz 
receiver. Several lines are detected suggesting 
the existence of multiple gaseous components to- 
wards the direction of CAHA J23056-^6016. The 
strongest emission line is seen at a velocity of rest 
of -49.3 km/s and a slightly weaker line at -51.3 
km/s (Table fl]) . The velocity of the second peak 
compares very well to the velocity we find for C^^O 
emission in our SMA data (see below). Thus, 
we believe that CAHA J23056-I-6016 is associated 
with emission at this velocity of rest and we assign 
a value of -51.3±0.2 km/s to the object. Using 



the galactic rotation model of Brand ( 1986p one 
finds a corresponding kinematic distance to the 
source of ~5.1 kpc. The detected velocity is fairly 
consistent with those of the nearby HII regions 
mentioned in section 3.1.1: For IC 1470 and IRAS 
23033+5951, iBronfman et al.|(|1996|) found veloci- 



ties of -53.1 and -52.1 km/s, respectively. Wouter- 
loot et al. ( 1989|) found for the later region a ve- 



locity of -51.8 km/s. Even if both HII regions lie 
-3.6' and -7.7' away from CAHA J23Q56-h6016 
one can assume that all belong to the same molec- 
ular cloud complex. In particular, Russeil et al 



(2007) concluded that their star-forming complex 
named "110. 1-1-00" (based on its galactic coordi- 
nates) has a velocity of -52.0 km/s and is com- 
prised of several HII regions including IC 1470 and 
also BFS 17, the third HII region close to CAHA 
J23056-I-6016 as mentioned in section 3.1.1. How- 
ever, these authors also revise the distance to indi- 
vidual HII regions based on the spectral type and 
magnitudes of the exciting sources of the regions. 
For IC 1470 and BFS 17, they derive distances 
of 2.87±0.75 kpc and 3.66±1.7 kpc, respectively, 
which is significantly less than the 5.1 kpc derived 
above from the galactic rotation model. Further- 
more, for IRAS 23033-^5951 |Reid fc Matthews] 

(2008) give a distance of 3.5 kpc which they as- 



sume for the whole Cepheus star-forming complex. 
To be in agreement with the most recent estimates 
for the other regions, we adopt a distance of 3.5 
kpc for CAHA J23056-t-6016. 

3.3.2. Disk mass estimate from, dust continuum, 

Additional mm-observations were carried out 
with MAMBOl at the IRAM 30 m telescope and 
finally with the SMA. The MAMBOl observations 
at 1.2 mm continuum yielded a flux of 12.9±1.0 
mJy with a beam FWHM of -10.5". Hence, the 
beam covered not only the dusty disk in the cen- 
ter but also the nebula. During interferometric 
follow-up observations with the SMA with a syn- 
thesized beam of —4" we measured a flux of 11.1 
mjy/beam at 1.3 mm. The almost identical flux 
levels of the single dish and the interferometric ob- 
servations show, that by far the bulk of the dust 
emission is coming from the central regions of the 
object, i.e. the disk. Since we did not resolve 

^^The model assumes a distance from the sun to the galactic 
center of 8.5 kpc and a velocity of the sun of 220 km/s. 



CAHA J23056+6016 with our SMA observations 
this emission comes from within the inner 4". Fig- 
ure |6] shows an overlay of the SMA mm-data on 
the iJ-band Subaru image. The peak of the mm 
dust emission coincides very well with the IRAC 
5.8 /^m flux and the dark lane seen in the NIR im- 
ages (Figure [2]). Hence, the mm data confirm that 
the dark lane is indeed a dusty disk. Since the 
disk is not resolved in our SMA data, we can only 
give an upper limit on its size based on the mm 
observations. For a distance of 3.5 kpc we find the 
radius of the dust disk to be < 7000 AU. 

Assuming on optically thin configuration, we 
obtain for the dust mass 



M, 



Fi. 



dust 



3 mm 



/^1.3? 



■B 



1.3 mm \-^ dust 



{Ti 



(2) 



where F is the integrated flux, d is the distance, 
K1.3 mm is the dust opacity in cm^/g at 1.3 mm and 
Bi.3mm(Tdust) is the Planck function at 1.3 mm 
for an assumed dust temperature Tdust- Since we 
do not know the temperature of the dust disk, nor 
the correct opacities for the grains, we estimate 
the total disk mass for a set of different param- 
eters in Table [2J Here, we always assume a flux 
of 11.1 mJy/beam and a distance of 3.5 kpc. The 
two values we use for the opacity correspond to 



that for dust grains with ice layers ( Ossenkopf k. 
Henning|1994 ) and that assumed for circumstellar 
disks around low-stars, where in general some dust 



processing is assumed to have happend (Andrews 
fc Williams||2005 ) . Furthermore, to derive the to- 
tal disk mass (gas and dust) from the dust mass, 
we multiply the values computed with equation (2) 
with the classical value for the gas-to-dust ratio of 
Mgas/Mdust = 100, although a factor of «190 has 
been derived as median gas-to-dust ratio of nearby 
galaxies by Draine et al. (2007). The estimated 



range of disk masses reaches from Mdisk = 0.4 M© 
to Mdisk — 2.7 Mq, reflecting the uncertainties for 
mass estimates if dust opacities and temperatures 
are not well constrained. 

3.3.3. Molecular outflow and hints for gas rota- 
tion 

In addition to the mm continuum data, C0(2- 
1) molecular line observations were carried out at 
the SMA simultaneously. In Figure |6] the green 
contours show the integrated C0(2-l) emission 



between -58 and -54 km/s. This blue-shifted elon- 
gated emission overlays nicely with the northern 
part of the nebula suggesting that the latter one 
is a cavity created by the molecular outflow. The 
fact that the line emission is blue-shifted sup- 
ports the hypothesis that the northern part of 
CAHA J23056-H6016 is shghtly tilted towards the 
observer. A red counterpart to this blue-shifted 
emission is not detected in the SMA data. This 
supports indirectly the hypothesis that the addi- 
tional CO emission lines seen at lower velocities in 
the SMT spectrum in Figure [5] are caused by fore- 
ground material and are not directly related to 
CAHA J23056-I-6016 because the emission is fil- 
tered out in the interferometric observations with 
the SMA. 

C^^ 0(2-1) emission is overplotted in blue and 
red contours in Figure [6) Although only detected 
in two adjacent velocity channels there is a clear 
spatial offset between the two peaks roughly in 
east-west direction and thus parallel to the dusty 
disk and perpendicular to the outflow mentioned 
above. These data, thus, seem to probe the 
gaseous counterpart to the dusty disk and the shift 
of the peak in velocity space can be interpreted as 
sign of rotation within the gas. It is interesting 
to note, that the gaseous disk component appears 
to have a larger physical extent than the dusty 
component probed in the continuum observations. 
However, there are certainly limitations to the cur- 
rent data set and the interpretation: While the red 
component of the emission appears elongated as 
expected for a rotating disk, the blue component 
is in addition much more extended in north-south 
direction. Furthermore the peak emission in the 
red and blue component is not perfectly aligned 
with the dark lane seen in the NIR images. Still, 
the C^^0(2-l) data support the picture of an ob- 
ject that is surrounded by a gaseous, rotating disk 
even if observations with higher sensitivity and 
spectral/spatial resolution are certainly eligible. 

3.3.4- Disk mass estimate from C^^ em,ission 

The C^^ 0(2-1) emission can also be used to 
estimate the related H2 column density and from 
this the related gas mass, if the distance to the 



object is known. Wilson et al. (2009) derive the 



H2 column density via 

iV(H2) = 2.65 • 10^^ Tdv^ 2.65 • 10^^ • / [cm^^] 

(3) 
where / is the integrated emission in K-km/s. In 
the Rayleigh- Jeans approximation we find / = 
6.86 K-km/s which translates into 



7V(H2) w 1.82 • 10^2 cm" 



(4) 



The mass of the gas can then be computed using 



Mgas = A^(H2) -A-mn-fi [g] 



(5) 



with A being the size of the emitting region, nin 
being the atomic mass unit, and /x is the mean 
molecular weight of the gas. To estimate the size 
of the emitting region we averaged the sizes of the 
blue and red contours in Figure |6] and found ^40 
arcsec^. Furthermore, we use mu = 1-66 • 10~^^ g 
and n = 2.7 to account for roughly 35% helium in 
the gas. With 1 Mq sa 1.99 • lO^^ g we derive 



M„ 



4.5 M, 







(6) 



This figure is at least a factor of 2 larger than the 
values we derived for the disk mass from the dust 



continuum data in section 3.3.2 However, as men- 
tioned above and seen in FigureM the gas emission 
is also more extended than the compact and un- 
resolved dust continuum emission. Together these 
estimates show that up to a few solar masses of 
gas and dust are confined in the immediate cir- 
cumstellar environment of CAHA J23056-I-6016. 

Although frequently done in the literature, we 
refrain from estimating the enclosed dynamical 
mass of the system by equating the gravitational 
force with the centrifugal force with the latter one 
being based on the observed rotation velocity and 
radius of the gaseous disk. This method implicitly 
assumes a Keplerian rotation of the disk which is 
no longer valid if the disk mass is not negligible 
compared to the central mass, which could be the 
case for CAHA J23056-f6016. In addition, there 
are examples for massive, rotating disk-like struc- 
tures that explicitly show non-Keplerian rotation 
(e.g., |Beutherfc"Walsh,2008j . 

3.3.5. HCO+/H^^CO^ excitation temperature, 
column densities and abundances 

With the IRAM 30 m telescope we searched 
for HCO^/H^^CO"'" signatures which are known 



to trace regions of higher density. The criti- 
cal density for the (1-0) transition at a temper- 



ature of 10 K is roughly 2 • 10 cm "^ (Evans 



1999[ ). Figure [t] shows that both lines, HCb+(l 
0) and H13CO+(1-0), were detected. At the 
HCO"'"(1-0) frequency, we find two lines separated 
by ca. 2.5 km/s. The stronger component is cen- 
tered roughly at the velocity found for C^^O in 
our SMA data. The much weaker H"CO+(1-0) 
line coincides also with the stronger HCO+(1-0) 
line at around -51.5 km/s. The nature of the weak 
HCO"'" emission at -49 km/s is unclear. Its rela- 
tion to CAHA J23056-t-6016 is unknown, and it 
might be a foreground contribution from material 
not directly associated with CAHA J23056+6016. 
Since we have not mapped the emission, we cannot 
further investigate this hypothesis at this point, 
but also in the CO (3-2) spectrum there are addi- 
tional emission components (see. Figure [5]). The 
-51.5 km/s emission on the other hand is most 
likely related to CAHA J23056-h6016 and demon- 
strates that dense molecular gas is associated with 
the object. The line strength ratio of the two iso- 
topologues (~ 6:1) indicates that the HCO+(1-0) 
emission is optically thick as typically the abun- 
dance ratio between ^^C and ^'^C is in the order 
of 60:1. Noteworthy are the narrow linewidths of 
0.7-1.0 km/s. Hence, the level of turbulence (po- 
tentially causing much broader lines) seems to be 
moderate at best. 



Following Purcell et al. (2006) we calculate an 



excitation temperature T^x for both lines and af- 
terwards the column density for the optically thin 
H^^CO"*" emission. Assuming that the HCO^ line 
is optically thick (1 — e""^ w 1), as suggested from 
the line ratios, we can calculate the corresponding 



Tex according to 



T — 

-^ ex — 



IWu 



In 1 






m ■ (7) 



Here, h is the Planck constant, k is the Boltzmann 
constant, and v^ is the rest frequency of the tran- 
sition (see, Table IT]). Furthermore, T^ is the beam 
corrected brightness temperature for HCO+ 



rri -^ antenna r-rj^i 

Vh ■ f 



(8) 



with Tantcnna being the antenna temperature (0.5 
K), 77b being the beam efficiency (0.85) and / being 



the assumed beam filling factor (0.05r^ Finally, 
we have 



^.(Tbg) 



hu„ 



1 

j/ii/u/fcTbg _ I 



(9) 



with Tbg as the temperature of the background 
radiation (2.73 K). Plugging in the corresponding 
numbers we find for the excitation temperature 



Te: 



13.8 K 



(10) 



and probably we can assume that the HCO+ and 
H-'^^CO"'" emission arises from the same gas and 
have hence the same excitation temperature. 

Knowing the excitation temperature, we can es- 
timate the optical depth of the optically thin line 
emission of H^'^CO"'" and also the column density 
of the molecule. The optical depth can be derived 
via 



Ttlu 



-In 



1 



Tex — Jy{Thg) 



(11) 



where T^ is now the beam corrected brightness 
temperature for H^^CO+ for which Tantenna = 
0.08 K. Beam efficiency and filling factor are the 
same as above. We find 



Tthv 



0.15 



(12) 



Assuming LTE conditions the total column 
density N can be estimated using 



7V = 7V, 



Q{Tex) 



1 



^-E^jkT^.. 



(13) 



where A'u is the column density of the upper state, 
g^ the degeneracy of the upper state, E'u the en- 
ergy of the upper state, and Q{Tex) the partition 
function. The upper-state column density can be 
computed via 



N^ = 



hc^A„ 



Tbdiy [cm"2] 



(14) 



with, again, k and h being the Boltzmann and the 
Planck constant, respectively, Aui being the Ein- 
stein coefficient of the transition in s~^, c being 
the speed of light, and ly^ and J_ Ti^dv being 



^*Based on the extent of the C^*0 emission in Figure^ we 
assume that the HCO"*" emission comes from within 6.5". 
With a beam size for the HCO"'' observations of roughly 
29" wc obtain a filling factor of 0.05. 
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the transition frequency and the measured hne in- 
tensity, respectively (see, Table [l]). The partition 
function QiT^^x) for linear rotators can be approx- 
imated by 



QiT,, 



hB 



(15) 



( Purcell et al.|20d6[ ) where B is the rotational con- 
stant i.e., 4 3377 MHz for H"CO+ ( |Schmid-Burgk 
et al.|2004 ). Plugging in the derived value for Tf 



and using g^ = 2J+1 — 3, A^i 
and Eu =4.16 K we derive 



3.8534-10" 



iV(H"CO+) « 7.0 • 10^° cm 



(16) 



for the H13C0+ column density yielding A^(HCO+) 
4.2 • 10^^ cm^^ assuming an abundance ratio of 
60:1 for i^C/^C. 

To get an estimate of the overall HCO+ and 
H^'^CO''" abundance we need to derive the H2 col- 
umn density. In section 3.3.4 we already esti- 
mated the H2 column density based on the in- 
terferometer data. For a direct comparison with 
the HC0+/H13C0+ IRAM 30m observations it 
is, however, more appropriate to derive the col- 
umn density also from single dish data. This 
can be done using the dust continuum emission 
we observed with the IRAM 30 m telescope as 
here the beams of the dust continuum and the 
HCO+/II^^CO+ observations are the most com- 
parable. It holds 

F ■ Rg/d r 



N{R2) = 



^ • B^(Tj_ust) ■ Hy-ma- ^ 



(17) 



see, e.g., Vasyunina et al.| 2009), where F is the 



dust continuum peak emission, Rg/d is the as- 
sumed gas-to-dust ratio, Vt is the solid angle (in 
steradians) covered by the beam, B^{Tdust) is the 
Planck formula for the assumed black body emis- 
sion of the dust at frequency v and temperature 
Tdust, i^u is the dust opacity, toh is the atomic 
mass unit, and /i is the mean molecular weight of 
the gas. For F we use 12.9 • 10~^^ erg/s/cm^/Hz 



(see section 3.3.2 ). Furthermore, we use Rgjd—^QQ 
and Ki,=1.0 cn?g~^. ttih and /i are the same as in 
section |3.3.4[ The beam solid angle can be com- 
puted via 



f^ = 27r('l-cos-][sr] (18) 

with uj being the beam size in radiarp] Although 

15 1 rad Ri 206265 arcsec 



the beam size for the IRAM30 dust continuum ob- 
servations was ~10.5", we have again to take into 
account the assumed size of the H^'^CO+ emission. 
Thus, to be consistent, we use an effective beam 
size of 6.5". Based on equation 17 we hence derive 
a hydrogen column density of 



iV(H2 



1.32 • 10^2 cm"2 



(19) 



yielding an abundance ratio for H^'^CO^ of 
[H13C0+] _ iV(Hi3C0+) 



X = 



[H2 



iV(H2 



5.3 • 10" 



(20) 



We will discuss these results below in section [531 
but note already that the H2 column density de- 
rived here is in general agreement with the value 
based on the C^^O SMA data (section [3A4). 



3.3.6. Detection of C2H 

In Figure Is] we show the C2ll(2-1) emission 
line towards CAHA J23056+6016 detected with 
IRAM30. The line is weak and very narrow (Ta- 
ble [1]) but clearly detected and the observed ve- 
locity agrees with the assumed velocity of rest of 
CAHA J23056-h6016 of -51.3 km/s. The possible 
origin and the interpretation of the emission will 
be further discussed in section [5.51 



4. Radiative transfer modeling 

The central illumination source in CAHA 
J23056-h6016 is obscured in the NIR and MIR 
by the circumstellar material, so that radiative 
transfer calculations are necessary to analyze the 
radiation that we receive from the dust particles 
near the source. 

Due to the scattering integral in the radiative 
transfer equation and the complex geometry of the 
dust configuration around the source, calculations 
of NIR scattered light images are time-consuming. 
Given the many free parameters that are neces- 
sary to describe the disk distribution, the radia- 
tion source, and the outflow, an automated fitting 
procedure and an extensive exploration of the pa- 
rameter space are therefore prohibitive. 

Moreover, the IRAC images, especially at 5.8 
//m, are expected to be influenced by PAH emis- 
sion which opens up an even wider parameter 
space. 
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For this paper, reporting on the detection of the 
new large, massive disk, we restrict the analysis to 
a simple model based on a 2D density distribution 
for the disk and the reflection nebula. 

From the images, it is obvious that a 2D mod- 
eling will not able to reproduce the warping of the 
disk and the H and K-h&nA flux maximum in the 
left part of the nebula over the dark lane. The 
model was chosen to estimate the disk and enve- 
lope parameters that are needed to reproduce the 
dark lane and the overall flux pattern in the re- 
flection nebula. 

The models make use of the constraints derived 
from the SMA image analysis, namely a total disk 
mass between 0.4 and 2.7 M0 (see Table 2), and a 
radial disk extent < 7000 AU (see section 3.3.2.). 
For the density distribution of the dust in the disk, 
we use a parameterized disk model with a radial 
power-law profile and an exponential disk atmo- 
sphere (e.g., Pascucci et al.|[2004 | of the form 

, 2" 



n(? 



"0 



''o 



exp 



\rh) 



(21) 



This dust number density is normalized by no, as- 
sumes a clear inner zone around the central ob- 
ject up to the inner radius tq, a truncation at the 
outer radius ri, the radial powerlaw exponent a, 
and the e-folding scale height ratio h = H/R {H is 
the z-value for which the density has dropped by 
1/e for a given radius). The disk may be inclined 
with respect to the edge-on view by the angle i, 
and rotated by S in the plane-of-sky relative to the 
north. 

To obtain a reasonable description of the 
boundaries of the reflection nebula, we normal- 
ized the H and K image equally and added them, 
rotated the resulting image to achieve a vertical 
symmetry axis, and added the image parts left and 
right from the axis. In this image, the border line 
of the nebula was fitted with a 4th-order polyno- 
mial. Figure [9] shows the resulting image as inlet 
and the contour plot with the fitting polynomial. 

For the density structure in the envelope, we 
have used several model types including a thin 
sheet at the nebula border and a density distri- 
bution filling the entire cavity. The finally used 
distribution has a maximum at the nebula border 
and a linear decrease toward the symmetry axis 
down to a constant density, as well as a linear de- 
crease along this axis. The distributions were at- 



tached at a height over the disk where the nebula 
border line height is equal to the disk scale height 
for the outer radius. The dust particles were as- 
sumed to be 0.1 ^m-sized spherical silicate grains 



using the optical properties described in Draine 



(2003) 



We have varied the luminosity of the central 
source to cover the range from young, interme- 
diate mass objects (Herbig Ae/Be stars, Herbij 
1960) to massive stars. The radiative transfer cal 



culations were performed using a ray-tracing code 
( [Steinacker et al.||2006[ ). 



4.1. Disk properties 

Figure [TO] shows the observed iiT-band image 
(left) and the synthetic image based on a disk den- 
sity model (right) with a dust number density of 
10^ m~'^ at the inner disk radius tq — 20 AU, a 
density powerlaw index of 1.8, a disk inclination of 
10°, a nebula tilting of 10° in the plane-of-sky, an 
outer radius of 4000 AU, and a disk scale height 
ration of 0.2. The outer radius in the model indi- 
cates that indeed the true disk size is smaller than 
inferred from the NIR images and that shadowing 
effects may play a role. The model was chosen by 
comparing the resulting image with the model im- 
ages of about 400 other models varying the disk 
parameters. We have restricted the modeling to 
the X-band image which shows the best signal-to- 
noise ratio. 

It reproduces the position, tilting, extent, and 
thickness of the dark lane. The asymmetric dark- 
ening in the right part of the disk is not repro- 
duced within the 2D model. The images are not 
sensitive to changes in the inner or outer disk ra- 
dius since the disk is seen almost edge-on and the 
outer disk parts are shadowed by the inner parts. 
The images are sensitive to the disk inclination 
since a larger inclination amplifies the scattered 
light flux from the upper disk part in disagreement 
with the observed flux pattern. Moreover, for in- 
clinations > 15°, the disk lane is stronger bend 
than observed. The model images also failed to 
reproduce the X-band image when the disk scale 
height deviated by more than 0.1 from the value 
0.2. There was no strong dependence on the radial 
spectral index. 

Since the outer disk parts do not modify the 
images strongly as long as the disk is kept al- 
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most edge-on, a variation in the disk mass with 
otherwise unchanged disk properties was possible 
within the mass range determined from the SMA 
data without changing the images in their overall 
appearance. 

The image shows stronger scattered light flux 
from the upper disk atmosphere compared to the 
flux from below the dark lane due to the inclina- 
tion of the disk. 

A clear difference between the observed and the 
synthetic image is the lower part of the reflection 
nebula which is much weaker than the upper part 
in the observed flux while the theoretical image 
shows comparable fluxes in the outer nebula parts. 
We have investigated if forward scattering or disk 
inclination could cause this effect. The scatter- 
ing phase function is too homogenous to reduce 
the flux in such a way, and a stronger disk incli- 
nation will produce too much flux from the up- 
per reflection nebula. Since this part shows also 
a stronger deviation from rotational symmetry in 
the observed image, we speculate that outer ma- 
terial may cause extinction patches which reduce 
the flux of the lower nebula. 

4.2. The central source 

A determination of the parameters of the cen- 
tral source is hindered by the ambiguity in the 
explanation of the scattered light flux from the 
nebula. The same flux value can be obtained when 
the stellar flux is decreased but the column density 
within the nebula is increased. Moreover, the dust 
opacities are poorly known, and will add another 
factor of 5 of uncertainty to the stellar proper- 
ties determined from the scattered light measure- 
ments. 

5. Discussion 

5.1. Disks around young stars - masses 
and sizes 

Putting the results from the previous sections 
in context with other young objects, the circum- 
stellar disk surrounding CAHA J23056-t-6016 ap- 
pears quite exceptional in two ways: (1) it is fairly 
large with an outer radius < 7000 AU, and (2) it 
is fairly massive with a mass in the order of up to 
a few M© 



(see, sections 



3.3.2 and 3.3.4) 



In par- 
ticular for the second point we emphasize that our 



SMA observations probe almost exclusively emis- 
sion from the disk and at least the mm-continuum 
data are not contaminated by additional flux from 
a remnant envelope. 

Circumstellar disks around young, low-mass 
stars (TTauri stars) typically have masses in the 
order of a few times 10~^-10~^ M0 and they have 
a maximum radial extent of a few hundred AU 



(e.g., Beckwith et al. 1990 Koerner & Sargent 



1995 Andrews & Williams 2005). In these sys- 



tems the disk-to-star mass ratio is roughly 0.5-1% 
although the ratio is certainly larger at earlier evo- 
lutionary stages when mass is continuously fed to 
the disk by a circumstellar envelope. For instance. 



Eisner & Carpenter ( 2006 ) found in a mm-survey 



in the Trapezium cluster the mean circumstellar 
disk mass to be in the order of the minimum-mass 
solar nebula (-0.005-0.01 Mq), but they found 
a few sources with disk masses between 0.13-0.39 
Mq (with at least a factor of 3 uncertainty) for 
which it seems likely that they are in an earlier 
evolutionary stage. 

For young Herbig Ae/Be stars the derived disk 
masses range from values also comparable to the 



minimum-mass solar nebula (e.g., Mannings & 
Sargent|[T997l |Thi et al"1[2004{ [XIo"nso-Albi et al 
2009[) to slightly higher m asses, i.e. -- 10"^ Mq 
(e.g., Henning et al.|1998 ). The sizes of these disks 
appear to be in the same range as those around 



lower mass stars or only slightly larger (e.g. Man- 



nings fc Sargent|1997[|Dutrey et al.|2007[|Schreyer 



et al. 2008 Alonso-Albi et al. 2009). There are 



at least two noteworthy exceptions, though: The 
first one is the disk around the pre-main sequence 
A-type star IRAS 18059-3211 where the gas disk 
extends out to 1800-3000 AU and where the disk 
mass is estimated to be in the order of 0.3 - 0.9 M© 



as derived from dust continuum data (Bujarrabal 



et al. 2008 2009 



Unfortunately the distance to 
this source is not extremely well constrained leav- 
ing these numbers with comparatively high uncer- 
tainties. In addition, an applied disk model pre- 
dicts a disk mass that is an order of magnitude 
lower than the values cited above. The second one, 
is the heavily discussed circumstellar disk in M17 



(e.g., I Chini et al.|2004HSako et al. 12005) [Steinacker 



et al.|2006[|Niirnberger et al.|2007[|Nielbock et al. 



2008 ) . Most recent results suggest that the central 



object is an intermediate mass star between 3-8 
Mq ([Nielbock et al.||2008[), but the disk was not 
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detected with the SMA and its mass and actual 
extent depend strongly on modeling assumption 
dSteinacker et al.||2006[ ). 

For high-mass stars (> 8 Mq), as mentioned 
already in the introduction, the existence of cir- 
cumstellar disks has not been fully established 
although there is growing observational evidence 
that they exist (e.g.JCesaroni et al.l'2007||Reid & 



Matthews|[2008l [Jiang et al.||2008i .Beuther et al 
20091). For example, iSchreyer et al'.l (|2006|) found 



evidence for a rotating ~1.5 Mq disk structure 
around the 8-10 M© object AFGL 490 based on 
interferometric observations of C^^O line emission. 
However, more often than not the main prob- 
lem is the lack of spatial resolution which would 
help to disentangle rotating "disk-like" structures, 
"pseudodisks" ( Galli fc Shu|1993 )and gaseous tori 
(possible remnants from the parental molecular 
cloud core) from real accretion disks surrounding 
and feeding material to the young massive star. 
But even with sufficient spatial resolution, e.g., 
using MIR interferometry around 10 /im where 
dusty disks are bright, so far the existence of large 
and/or massive circumstellar disks around mas- 
sive young stellar objects could not be proven 
(iLinz et al.||2009l Ide Wit et all 120091). Only re- 



cently Okamoto et al. ( 2009 ) found direct evidence 

10 M0, young star 
This radius 



for a flared disk around the 
HD200775 from direct MIR images, 
of the disk extends out to 1000 AU, but the mass, 
derived from unresolved emission at 350 GHz, is 
only a few 10~^ M©, with the exact value depend- 
ing on the assumed dust opacities and tempera- 
ture. Indirect evidence for accretion disks around 
young objects of similar mass has been found by 



Boley et al. (2009) and Davies et al. (2009) for S 



235 B and W33A, respectively. Resolved images 
of the latter two disks have not yet been obtained. 

This shows, that CAHA J23056-H6016 is a good 
candidate for being among the largest and most 
massive circumstellar disks detected and resolved 
on multi wavelengths observations. However, one 
should keep in mind that it has not yet been shown 
that the disk is indeed an accretion disk from 
which matter is transferred onto the central ob- 
ject while angular momentum is transported out- 
wards. The initial concept of "pseudodisks" ac- 
tually predicts gaseous, flattened structures with 
sizes at least comparable to what is observed here 



the current instrument sensitivities the dusty disk 
component of CAHA J23056-f 6016 appears to be 
more compact than the detected gaseous disk com- 
ponent (section 3.3.3, Figure [6]). It would be inter- 
esting to check observationally for signs of accre- 
tion coming from a compact inner accretion disk 
that is embedded in a somewhat larger flattened 
molecular structure. 

5.2. Mass of the central object 

Our current data does not allow us to constrain 
the mass of the central source empirically. Only 
indirectly, by comparing the derived disk mass to 
those of other systems as done above, we can spec- 
ulate that most likely the central source of CAHA 
J23056-I-6016 is at least an intermediate mass ob- 
ject with several solar masses. Trying to derive 
the spectral type of the central source from NIR 
spectroscopy of the scattered light in the reflection 
nebula could possibly be a means of getting more 
direct information of the object. 

5.3. CAHA J23056+6016 - an evolved ob- 
ject? 

Up to now it was assumed that CAHA J23056-I-6016 
is a young stellar object. But what makes us be- 
lieve that it is indeed young and not evolved? 
At first glance, the NIR morphology of CAHA 
J23056-I-6016 is not only compatible with a young 
object but does also resemble that of proto- 
planetary nebula (PPNe), i.e., post-AGB stars 
(e.g., see Meakin et al. 2003 for high-resolution, 
NIR HST images of OH 231.8-^42). In addi- 
tion, also for the intermediate mass, pre-main 
sequence object IRAS 18059-3211 (already men- 
tioned above) Ruiz et al. (1987) initially suggested 
that the object is rather an evolved source and not 
a pre-main sequence star. 

We note that despite the apparent similarities 
in the morphology compared with evolved stars, 
CAHA J23056+6016 is lacking some crucial prop- 
erties of older objects. (1) [Meakin et al. (2003) 



(Galli & Shu 19931. Interestingly, however, with 



found a very filamentary sub-structure in the bipo- 
lar cones surrounding the PPNe OH 231.8-1-42. 
This structure is caused by material that is ejected 
from the central source. This mass loss is, how- 
ever, neither constant nor homogenous in space 
or time and leading to knots and filaments in the 
bipolar cones. Such a sub-structure is not seen 
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in the nebula around CAHA J23056+6016 (Fig- 
ure [2|. The nebula around CAHA J23056+6016 
appears rather smooth and its brightness distribu- 
tion can be very well explained with an illuminated 
outflow cavity with small dust grains scattering 
the light coming from the central source. (2) The 
mass loss of evolved stars is typically accompanied 
by high- velocity molecular outflows and jets (e.g., 
Sanchez Contreras & Sahai 2004) and the observed 



outflow velocities of PPNe typically exceed several 
tens of km/s with respect to the system's veloc- 
ity of rest ( [Bujarrabal et ah 2001). But, again, 
we do not find any evidence for a coUimated jet 
or a high-velocity outflow component for CAHA 
J23056+6016, neither in the NIR images (where 
shock-induced hydrogen emission could show up 
in the K'-hand image) nor in the CO data. (3) Fi- 
nally, evolved objects also show expanding molec- 



ular tori lying in the disk plane (e.g., Hirano et al 



2004 Meakin et al. 2003). In the case of CAHA 



J23056-I-6016, however, the molecular gas that we 
see in our C^^O data is well confined to a small ve- 
locity regime making the interpretation of a rotat- 
ing, disk-like structure much more likely compared 
to an expanding torus. 

Based on these findings and the fact that we 
are confident that CAHA J23056+6016 is indeed 
a young and not an evolved object. 

5.4. An isolated massive circumstellar 
disk 

In one way CAHA J23056-I-6016 appears to be 
different compared to the majority of young in- 
termediate mass and higher mass objects. Most 
of the young, sources mentioned above happen to 
be in young clusters or star-forming regions (e.g., 
M17, HD200775, S 235 B, W33A, as well a s the 
lists of massive protostellar disk candidates in Reid| 
k Matthews! (|2008|) and [Cesaroni et all (|2007[)). 



For massive stars, Smith et al. (20091 (among oth- 
ers) claim that they can only form in cluster en- 
vironments and that it is the overall evolution of 
the cluster which leads to the formation of mas- 
sive objects in the center eventually. Examples for 
this can be found in Figurefll where one can nicely 
see how the young, more massive stars in the HII 
regions are surrounded by a noticeable group of 
additional objects. The sizes of these clusters are 
in very good ag-reement with the findings from 



young stellar cluster to be in the order of 0.5-1 
pc which corresponds to roughly 0.5-1' in Fig- 
ure [T] for a distance of 3.5 kpc. In this context it 
is even more striking that CAHA J23056-I-6016 is 
rather isolated while in other regions of the same 
molecular cloud complex star formation occurs in 
a cluster-mode. 

One interesting exception to the rule appears 
to be the intermediate mass, pre-main sequence 
star IRAS 18059-3211, which, as noted above, 
was initially classified as an evolved and not a 
young source. Indeed, the vicinity of CAHA 
J23056-h6016 looks very similar to that of IRAS 
18059-3211 as there is no direct evidence for a sur- 
rounding young cluster (Figure p] and Figure E| . 
Based on comparison to models for stellar evolu- 
tion of low-mass stars (jBaraffe et al. 



1998) 



estimate that in our deep NIR images we are able 
to detect objects with masses below 0.1 Mq for 
an assumed age of 1-5 Myr at a distance of 3.5 
kpc even if we assume a line-of-sight extinction 
of Av=l5 mag. However, we do not find direct 
evidence for a (tight) cluster of young sources of 
lower mass. Also, the reddish nebula south of the 
brightest star lying west of CAHA J23056-h6016 
in Figure [4] does not harbor any young embed- 
ded source: We did not detect this object in di- 
rect mm-continuum observations with the IRAM 
30 m telescope carried out in the same night as 
the observations of CAHA J23056-t-6016 (see sec- 
We rather assume that the nebula 



2.4.2) 



tion 

is reflective in nature and associated with the 
brighter star lying north of it. This star and 
the bright object east of CAHA J23056-K6016 in 
Figure |4] are the brightest NIR source in the di- 
rect vicinity of C AHA J23056-f601 6. Based on 
their NIR colors (Cutri et al. 20031) we estimate 



their spectral types to be K2-K4 (Ducati et al 



2001] ) which would make them foreground objects 
compared to CAHA J23056-h6016. Thus, even 
if we cannot exclude that maybe multiple low- 
mass sources are hidden behind the dusty disk 
and/or the nebula around CAHA J23056-h6016, 
we currently have no indications for a young clus- 
ter within 0.5 - 1 pc (the typical cluster size) of 
CAHA J23056-h6016. If this was the case, a solid 
determination of the central mass would be cru- 
cial. As mentioned above, in one theory more 
massive objects are thought to form mostly (if not 



Lada & Lada ( 2003 1 who found the typical size of 



solely) in company with other objects (Bonnell & 
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Bate 2006 Smith et al. 2009). The detection of 



CAHA J23056+6016, if it was a massive object, 
would talk very much in favor for a star forma- 
tion process for massive stars that is very similar 



to that for low-mass object (McKee & Tan 2003 
p<!rumholz et al. 2009): Single objects can form 



from the collapse of a molecular cloud core, and 
to reach its final mass an object accretes material 
from a surrounding circumstellar disk. 

5.5. Chemical abundances 

5.5.1. H^^CO+ andHCO+ 



In section 13.3.51 we derived an H^'^CO+ abun- 
dance of roughly 5.8 • 10~^^ which yields an HCO+ 
abundance of 3.5 • 10~^° assuming a ^^CO/^^CO 
ratio of 60:1. These values can be compared to 
those found in young cores and envelopes and also 
circumstellar disks. 

For two massive embedded protostars in the 



IRAS 18151-1208 region [Marseille et al.| ( |2008[ ) 
found an Hi3C0+ abundance of 3 - 8 • iF^. 



Similarly, Zinchenko et al. (2009) found for their 



census of regions of massive star formation val- 
ues between 5 • 10~^^ and 3 • 10^^°. Interest- 
ingly, for the envelopes of low-mass protostars 



van Kempen et al. (2009) found values that are 



roughly two to three orders of magnitudes higher 
(2 - 8 • 10-8). |Thi et al.| (|2004[), however, found 



for the solar-mass protostar IRAS 16293-2422 an 
H^'^CO''" abundance of '-^ 2 • 10"^^ which is similar 
to the values mentioned above for high-mass pro- 
tostars. The same authors provide also H^'^CO+ 
abundances for circumstellar disks around low- 
mass stars. Here, the values appear to be even 
lower, although a large dispersion exists. For 3 
low-mass disks, values between ^ 3 • 10^^'^ and 
~ 3 • 10"^^ were found. In addition, for 2 disks 
around Herbig Ae/Be stars Thi et al. (2004) give 
HCO+ abundances of 1 • 10^^'"^ 
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Given the large spread in these values, it seems 
difficult to derive any firm conclusion about the 
scientific meaning of the abundances we find for 
CAHA J23056-f 6016. In addition, the production 
and depletion of the molecules depend strongly 
on the environment (iThi et al.||2004|). HCO+ is 



mainly produced by HJ -I- CO — > HCO -f H2 in 
a gas phase reaction and here the formation is, 
among other parameters, dependent on the sur- 
rounding ionization degree. For the depletion. 



there are in general different mechanisms possible, 
e.g., dissociative recombination (HCO"*" -I- e~ — >■ 
CO -F H), reaction with water (HCO+ 4- H2O -> 
HaO"*" -I- CO) or charge transfer with low ioniza- 
tion potential elements such as Na or Mg ( Nikolic 



2007 ) . Again, these reactions depend on the envi- 



ronmental conditions. 

Despite all these uncertainties, we can note at 
least two things: First, the abundances derived 
for CAHA J23056-t-6016 lie in the same range as 
those found in circumstellar disks around young, 
low-mass or intermediate mass stars, but seem 
lower than the values found in (high-mass) pro- 
tostars. Second, using the derived abundances as 
lower limit for the ionization fraction in the disk, 
we find that the value is high enough to cause mag- 
netorotational instabilities in the disk which then 
drives turbulence and mixing (e.g., [Balbus^^2009| ) . 
Hence, we can expect a certain amount of disk 
dynamics to occur and possibly also accretion if 
the central source features a magnetic field that 
suffices to reach and tap the ionized parts of the 
disk. 

5.5.2. CoH 



As seen in section 3.3.6 we detected C2H emis- 



sion towards CAHA J23056-I-6016. The detection 
of this molecule and its interpretation in the as- 
trophysical context of star formation is subject of 
very recent publications in the literature. 



Beuther et al. (2008) found that C2H was 



present in a larger sample of massive star-forming 
regions covering various evolutionary stages in 
the star formation process. They suggest that the 
molecule could potentially be used to trace the ini- 
tial conditions of massive star formation, as in the 
earliest phases its abundance remains high due to 
constant replenishment of elemental carbon from 
CO being dissociated by the interstellar UV pho- 
tons. In later stages, the molecule is transformed 
into other species in the dense protostellar core 
regions and might only remain traceable in the 
outer core regions where it is excited by external 
UV radiation. 

In circumstellar disks, there is only evidence of 
C2H in three cases (DM Tau, LkCal5, and MWC 
480; Henning et al. 2010, submitted). However, 
the detection of the molecule in the disk around 
the Herbig Ae star MWC 480 was very weak. 
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From their observations and modeling the chem- 
ical evolution of the disks the authors concluded 
that C2H is indeed a sensitive tracer of the UV 
radiation field and that the stellar UV luminosity 
affects the evolution of C2H in the disk globally, 
while X-ray radiation appears to be more impor- 
tant in the inner disk regions close to the star. In 
this way, the lack of C2H in the Herbig Ae disk 
compared to the TTauri disks is explained with 
the higher UV flux from the central source. 

In this context, we suggest that the C2II we 
detect towards CAHA J23056-K6016 might be an- 
other hint for the young age of the source. How- 
ever, the current observations lack the required 
spatial resolution to conclude where the emission 
is coming from. Apart from emission from the 
disk itself it seems possible that the C2H flux is 
coming from the outer and thus externally illumi- 
nated parts of the disk / nebula system and not 
necessarily from the inner regions. In this case the 
nearby HII regions (Figure IT]) could be the source 
of the required UV photons to excite the emission. 

6. Summary and conclusions 

In this paper we presented a new large and mas- 
sive circumstellar disk (named CAHA J23056-t-6016) 
presumably surrounding a young stellar object. 
Our findings can be summarized as follows: 

• The disk is seen almost edge-on and was de- 
tected in high resolution NIR images where 
it caused a dark lane intersecting a bipolar 
nebula that is most likely illuminated by a 
young stellar source. The dark lane stretches 
over ~ 3" which, for an assumed distance of 
3.5 kpc, translates into ~10500 AU, but due 
to shadowing effects the true physical disk 
size is probably smaller (see also last point). 

• CAHA J23056-H6016 appears to be one of 
the largest circumstellar disks known today 
but at the moment we have no empirical ev- 
idence whether it is an actual accretion disk 
or a " pseudodisk" resulting from the collapse 
of a moderatly magnetized rotating molecu- 
lar cloud core or a flattened remnant enve- 
lope. 

• The young stellar source residing behind the 
disk is not detected in the NIR images as the 
optical depth through the disk is too high. 



• In Spitzer IRAC images the object is also 
detected and at 5.8 /xm one starts probing 
directly through the disk and detects flux 
mainly coming from the central disk regions. 

• This finding is supported by the fact that the 
peak of the emission seen in different Spitzer 
IRAC filters is shifting with wavelength and 
that at the longest available IRAC wave- 
length (5.8 /im) the emission peak coincides 
with the center of the dark lane seen in the 
NIR images. 

• In interferometric mm dust continuum ob- 
servations with the SMA we detected but 
not resolved CAHA J23056-I-6016 at the 11.1 
mjy/beam level. This sets the upper limit 
for the disk radius to < 7000 AU. 

• The mm dust emission peaks again nicely 
on the dark lane seen in the NIR images and 
coincides perfectly with the IRAC 5.8 /im 
emission. From the unresolved continuum 
data we estimate the disk mass to be in the 
order of 1 Mq making it one of the rare cases 
where a massive disk is resolved in NIR im- 
ages. 

• We found a molecular outflow in C0(2-l) 
coinciding perfectly with the northern part 
of the bipolar nebula seen in the NIR im- 
ages. In addition, we found indications for a 
gaseous disk component in C^^O data, with 
the emission being perpendicular to the out- 
flow and showing signs for rotation. The 
gaseous disk is even larger than that derived 
from the NIR images or the dust continuum 
observations and yields a mass of ^-^4.5 Mq. 

• We also detected Hi3C0+ and HCO+ line 
emission from CAHA J230564-6016 which 
arises from regions of higher density. The 
derived abundances are comparable to those 
found in other circumstellar disks but lower 
compared to protostellar source. 

• The detection of C2H supports a young age 
for the object. We suggest that the emis- 
sion could either arise from the disk itself but 
could also be emitted from externally illumi- 
nated regions of the disk / envelope system 
that are subject to high UV radiation flelds 
from nearby HII regions. 



17 



• The current data does not allow us to con- 
strain the mass of the central object directly 
but, based of the derived disk mass, we spec- 
ulate that it could have up to several solar 
masses as well. 

• Interestingly, we find no signs of a young 
cluster around CAHA J23056-I-6016 within 
the typical cluster size of 0.5-1 pc. If the 
central source were rather massive as well 
then we had found a rare case of isolated 
massive star formation. 

• Using a simple 2D dust density profile for 
the disk and the reflection nebula, the radia- 
tive transfer modeling of the i^-band image 
could constrain the inclination of the disk to 
be 5-15° from an edge-on view. The scale 
height ratio of the disk ranges from 0.1 to 0.3 
to be consistent with the extent of the dark 
lane. The radial power-law index was not as 
well constrained as the inner and outer ra- 
dius of the disk, where the latter parameter 
is found to be -4000 AU. 

The disk of CAHA J23056-I-6016 is extremely large 
and very massive and the fact that we have spa- 
tially resolved information at a variety of wave- 
lengths makes this object an ideal laboratory to 
study the properties of a disk surrounding a young 
presumably intermediate mass or even massive 
star. Our finding argues for a star formation pro- 
cess that is functioning over a wide mass range, 



ranging from sub-stellar objects (e.g., Zapatero 
Osorio et al. 2007 Quanz et al.| 2010 1 to, as it 
appears, at least several solar masses. In addi- 
tion, as higher mass stars seem to process their cir- 
cumstellar disks faster than low-mass objects (e.g.. 



Alonso-Albi et al. 2009 ), CAHA J23056-h6016 pro- 



vides the opportunity to study an early stage of 
disk evolution where the disk is still massive. 

For future observations of CAHA J23056-)-6016 
it might be interesting to aim at an even higher 
spatial resolution of the dust continuum emission 
to eventually resolve the disk. Filling in the SED 
in sub-mm and MIR wavelengths could in addi- 
tion help to even better characterize the disk ge- 
ometry and dust properties. In parallel, observa- 
tions of the gaseous disk component with higher 
spectral resolution will shed additional light on the 
detected rotation signatures. Finally, it might be 



worth to obtain a NIR spectrum of the reflection 
nebula. This might allow us to search for spectral 
features that constrain the spectral type of the 
central source and possible bear signs of accretion 
processes. 
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Table 2 

Estimated disk masses in units of solar masses for different values for the dust opacity k 

and different dust temperatures. a gas-to-dust ratio of 100 was assumed. 



K at 1.3 mill 


Td„,t=20 K 


Tdust^^O K 


Ref. for opacity 




2.7 Mq 
1.2 Mq 


0.9 M© 

0.4 Mq 




1.0 cmVg: 


Ossenkopf & Henning' ('1994|) 


2.3 cmVg: 


Andrews & Williams^ (2005) 
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Fig. 1. — NIR color composite consisting of J 
(blue), H (green) and Kg (red) filter images ob- 
tained at Calar Alto with Omega2000. The bipo- 
lar nebula intersected by a dark lane at the image 
center is CAHA J23056-f6016 (north is up, east 
to the left). 
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J filter 



H filter 



Ks filter 






Arc Seconds 



Arc Seconds 



Arc Seconds 



Center: R.A. 23 05 37,14 Dec +60 15 48,4 Center: R,A, 23 05 37,46 Dec +60 15 48.5 Center: R.A, 23 05 37,46 Dec +50 15 48,5 



Fig. 2. — Contours of the NIR images: J filter 
(left), H filter (middle), and K' filter (right). In 
the J band the S/N is lower than in the other 
filters and the 5 contour levels start at 5-cr and 
go up to I3-cr in equidistant steps. The lowest 
level corresponds to ~ 2.75 • 10"'* mJy/pixel. In 
the H and K' filter images the S/N is higher and 
the flux density levels of the contours are identical 
for better comparison. Here, the lowest contour 
equals ~ 2.0 • 10""* mJy/pixel corresponding to 
^1-a in H and ^Q-a in K' . The contours increase 
by 2.0 • 10^"* mJy /pixel reaching thus ^^ 6.0 • 10~^ 
mJy/pixel in the third level (maximum level in 
the H image) and ^ 1.0 • 10^'^ mJy/pixel in the 
fifth level (maximum level in the K' image). The 
triangle, the square and the cross in the H and K' 
image denote the position of the flux peak of the 
IRAC 3.6 /im, 5.8 ^m and SMA 1.3 mm emission, 
respectively. The point-like structure seen in H 
and K' ^ 1" north of the image center has very 
red colors and is presumably a background source 
shining through the nebula. 
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IRAC 1 



RAC 2 





IRAC 3 




3 2 1 0-1-2 

Arc Seconds 
R.A. 23 05 37.25 Dec 



3-4 3 2 10-1-2-3-4 3210-1-2 -3 -4 

Arc Seconds Arc Seconds 

f60 15 48.5 Center: R.A. 23 05 37.25 Dec +60 15 48.5 Center: R.A. 23 05 37.25 Dec +60 15 48.5 



Fig. 3. — Contours of the IRAC images: 3.6 /im 
(left), 4.5 /xm (middle), and 5.8 jim. (right). The 
lowest contour of the 3.6 fim image denotes a 6- 
a detection corresponding to ^0.0144 mJy/pixel. 
The contour levels increase in 2-a steps yielding a 
12-a level for the highest contour. For the 4.5 /xm 
image the contours start also at the 6-a level corre- 
sponding to ^0.0168 mJy/pixel but then increase 
in 3-CT steps with a maximum of 18-cr. Finally, 
in the 5.8 /.tm image the emission rather strong 
possibly enhanced by diffuse PAH emission. The 
contours denote 33-a at the lowest level equaling 
~0.099 mJy/pixel and increase in 3-a steps to a 
maximum level of 45-cr. 
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CENTRE: R.A. 23^05"'37?56 DEC +60°15'44'.7 (J2000) 



Fig. 4. — Color composite of -ff (blue), K' (green) 
and IRAC 4.5 /im (red). While the bipolar nebula 
is only seen in scattered light in the H and K' 
band, the flux in the IRAC filter nicely traces the 
disk intersecting the nebula. 
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Fig. 5. — Main beam brightness temperature of 
CO(3-2) transition line emission towards CAHA 
J23056+6016 measured with the 345 GHz receiver 
at the SMT. The parameter for gaussian fits to the 
different fines are given in Table IT] 




] Grey: Subaru H-^and 
Black contours: SMA 1.3mm continuurrt,; 
Red/blue: SMA c"o(2-1) 
Green: SMA btue-shifled C0(2-l) 
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Fig. 6. — Overlay of mm continuum and molec- 
ular line observations on the H-Sitei image. 
Black contours are 1.3 mm continuum emission, 
green contours blue-sfiifted CO (2-1), red/blue are 
C^*0(2-l). The continuum contours start at 3- 
a and increase in 1-a steps corresponding to 1.5 
mjy/beam. The green contours show the in- 
tegrated C0(2-l) emission between -58 and -54 
km/s. They start at 3-cr and increase in 1-a steps 
corresponding to 66 mJy/beam. The blue and 
red contours are centered on -51.6 and -51.0 km/s 
in velocity space, respectively, and each velocity 
channel has a width of 0.6 km/s. The lowest 
contour shows again the 3-cr level and the con- 
tours increase in 1-cr steps corresponding to 144 
mJy/beam. 
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Fig. 7. — Main beam brightness temperature of 
the HCO+(1-0) and Hi3CO+(l-0) emission lines 
towards CAHA J23056-K6016. The fit resufis 
(with set weight sigma) are summarized in Tablell] 



Tmb[K] 




-50 

Velocity (km/s) 



Fig. 8. — Main beam brightness temperature 
of the C2H(l-0) emission line towards CAHA 
J23056-I-6016. The fit parameter are summarized 
in Table [H 
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X [1000 AU] 



Fig. 9. — Upper reflection nebula borders derived 
from a convolution of the H and K images assum- 
ing symmetry. The inlet shows the convolved im- 
age, the contour plot contains the fit of the border 
by a polynomial as thick line. 




Fig. 10. — K-h&nd image (left) and scattered light 
model image (right) of the reflection nebula pro- 
duced by a radiative transfer calculation based on 
a 2D disk plus nebula density model. 
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